Transport of Entanglement Through a Heisenberg-XY Spin Chain 
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The entanglement dynamics of spin chains is investigated using Heisenberg-XY spin Hamiltonian 
dynamics. The various measures of two-qubit entanglement are calculated analytically in the time- 
evolved state starting from initial states with no entanglement and exactly one pair of maximally- 
entangled qubits. The localizable entanglement between a pair of qubits at the end of chain captures 
the essential features of entanglement transport across the chain, and it displays the difference 
between an initial state with no entanglement and an initial state with one pair of maximally- 
entangled qubits. 
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Quantum entanglement in spin systems is an 
extensively-studied field in recent years, in the advent 
of growing realization that entanglement can be a re- 
source for quantum information processingii^. Within 
this general field, entanglement of spin-i degrees of free- 
dom, qubits, has been in focus for an obvious reason 
of their paramount importance for quantum computers, 
not to mention their well-known applicability in various 
condensed-matter systems, optics and other branches of 
physics. For a pure state of many qubits, quantum en- 
tanglement is a measure of how a subsystem is correlated 
to the rest of the system, which can be quantified by the 
von Neumann entropy of the reduced density matrix of 
the subsystem. 

For a system of a large number of qubits, how different 
pairs of qubits are entangled in a pure state cannot easily 
be specified, even with known diagonal and off-diagonal 
correlation functions. There are many pairwise entangle- 
ment measures of pure quantum states of many qubits, as 
we shall discuss below. Essentially, any two-qubit entan- 
glement measure depends on the reduced density matrix, 
Pij of the marked pair of qubits i and j. The reduced 
density matrix is obtained from the original many-qubit 
pure state through a partial trace over the rest of 
the qubits, pij = Tr\^ij]^\'ip){ip\. Using the two-qubit basis 
|00), |01), |10), the reduced density matrix has the 
form 



^tj Hij 

^tj Cy 



(1) 



\ Gj^j I^j J^j Dij J 



The two-qubit reduced density matrix above represents 
a mixed state in general, though we started with a pure 
many-qubit state. We can calculate von Neumann en- 
tropy from the eigenvalues of pij, which quantifies how 
these two sites are entangled with the rest of the system. 
However, a measure of entanglement between these two 
sites cannot be straightforwardly quantified, because of 
the mixed-state structure. In the above each of the ma- 
trix elements can be expressed as an expectation value 
of a two-qubit operator in the initial state \ip). Using 
the Pauli operators for a qubit, cr+|0) = |1),(t+|1) = 



0,cr^|0) = -|0),cr^|l) = for instance, 
if,, = (01|p,,|10) = (V^|a+a,-h 



(2) 



and similarly the other matrix elements are expectation 
values of appropriate operators. A good measure of im- 
purity of the two-qubit state is Mij = 1 — Trpf^ , which 
is zero if represents a pure state of the given pair of 
qubits. In general the reduced density matrix represents 
a mixed state, giving a nonzero value for My . The von 
Neumann entropy calculated as —Trpij log pij quantifies 
entanglement of this pair with the rest of the qubits. 
When Mij = 0, the reduced density matrix represents a 
pure state for the pair of qubits and the von Neumann 
entropy would be zero. For Mij ^ 0, the von Neumann 
entropy will be nonzero, indicating an entanglement of 
the pair with the rest of the system, and possibly a mu- 
tual entanglement among the pair of qubits. The pairwise 
concurrence^ is a measure that quantifies the mutual en- 
tanglement of the marked pair of qubits. It carries the 
information of entanglement sharing between site i and 
any other site. If the concurrence between i and j is one, 
would imply neither of these two sites is entangled with 
the rest of the system, and hence a pure-state assignment 
can be done to these sites. If the concurrence less than 
one, that would imply site i is sharing entanglement with 
site j and the rest of the system too. The concurrence 
measure^ is given as 



Cij = max (0, Xi 
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In the above Xi are the eigenvalues in decreasing order 
of the matrix PijPij, where pij is the time- reversed ma- 



trix, /5,y 



y * 



<7:. The above concurrence 



measure captures the essential information of mutual en- 
tanglement of the marked pair of qubits. However, it is 
not clear if a given value of the concurrence for a pair 
of qubits would imply that one can separate the pair of 
qubits in an entangled state, and use the entanglement 
as a resource operationally for quantum information pro- 
cessing. 

The localizable entanglement (LE) between two qubits 
of a many-qubit state is the maximal amount of entan- 
glement system that can be concentrated on the given 
pair of qubits, on an average, by doing a measurement 
on the rest of the qubits^. 
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The LE between two qubits can be isolated, by an ap- 
propriate measurement on the rest of the qubits, and thus 
can be used as resource. Also, this has a direct bearing on 
the correlations present between the pair of qubit in the 
parent many-qubit state. The calculation of LE could be 
very tedious, as one has to search through every possible 
measurement basis for the rest of the qubits. However, 
there are useful upper and lower bounds on LE, that are 
calculable from the diagonal correlation functions alone^. 
Let us denote 



(4) 



and similarly Qfj , Q^j in terms of yy and zz correlation 
functions, and 



(l±(afa|)f-(«)±(a|)f 



(5) 



In terms of the above correlation functions, the bounds 
on LE are 



max {\Q^^\,\Ql\ 



\QW<LE, 




(6) 



Now, each of the matrix elements of the reduced den- 
sity matrix pij involves correlation functions. After some 
manipulations using the properties of the Pauli matri- 
ces, and their expectation values in the initial state, we 
have = 2Re(iJ + G) - 4:Re{E + J)Re(/ -I- F), Qy = 
2Re(i/-G)-4Im(£; + J)Im(/ + F),Q^ = 4:{AD-BC), 



and 



16AD, s- = 16BC. Here the ij subscripts 



have been suppressed for clarity. 

The entanglement properties depend only on the struc- 
ture of entanglement of a many-body state without ref- 
erence to a Hamiltonian. However, here we would like 
to investigate the entanglement dynamics of a simple in- 
teracting spin systems on a lattice. Consider a spin-i 
spin system, a nearest-neighbour anisotropic Heisenberg 
model in one dimension with a Hamiltonian 



(7) 



where Sf is spin-raising operator at site K is the 
exchange interaction strength, and A a measure of 
anisotropy. A constant, Eq the energy of a state with 
all down spins or all up spins, has been added on for 
convenience. Working in a diagonal basis of Sf for every 
site i, there are two states per site, viz., an up-spin state 
(denoted by |1) and a down-spin state denoted by |0)). 
A many-spin state can be characterized by the number 
of down spins, as the total 5^ is a good quantum num- 
ber. We confine ourselves here to spin systems in one 
dimension, as eigenstates are known in this case through 
Bcthe-Ansatz, and are well studied. 

The concurrence measure of entanglement in the eigen- 
states of the above model have been investigated^iSii, and 
the entanglement dynamics have been discussed^iSiifi. 
The localizable entanglement has been studied for the 
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FIG. 1: The entanglement measures for the last pair of qubits 
in an unentangled initial state \ipu) are plotted vs. time, for 
a finite open chain of A'' = 20 qubits. The upper bound on 
LE coincides with the concurrence measure in this sector with 
exactly on down spin. In the second figure the entanglement 
for nearest neighbour qubits 1,1 + 1 are plotted vs / for two 
different times Kt/h = 50, 150. 



ground stateW . In this study we will focus on the dynam- 
ics of the various entanglement measures for the above 
model. Our strategy is to start from a given initial state 
(either with no entangled pairs of qubits or with exactly 
one pair of qubits maximally entangled and the rest in 
an unentangled state), and study its entanglement prop- 
erties using the measures discussed above, through the 
time evolution using Heisenberg dynamics. The initial 
state |'0(O)), let us say with exactly one entangled pair 
of qubits at one end of the chain, will evolve into a state 
li'it)) at a later time t, with a complicated distribution 
of entanglement, with dynamically changing correlation 
functions. We shall investigate how entanglement trans- 
port takes place across the spin chain, by studying the 
entanglement of a pair of qubits at the other end of the 
chain. 

The dynamics of Heisenberg Hamiltonian above con- 
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serves the number of down spins in a given state. This 
imphes that, for an initial state with a definite number for 
down spins, at future times the reduced density matrix 
for a pair of qubits can have nonzero diagonal elements 
and a nonzero off-diagonal element Htj in Eq. (1), and 
we have 



(8) 



which simplifies the expressions for the entanglement 

measures substantially. We have for a many-qubit state 
with a definite number of down spins (number of qubits 
in state |0)), 



Mij = 
C - = 
LEij > 
LEij < 



2max(0,|i/y | - y^A~D~), 
max {4:\Aij Dij - Bij Qj \ , 2ReHij ) . 



(9) 
(10) 

(11) 
(12) 



A further simplification occurs for a state with exactly 
one down spin. Here, Aij the first diagonal element of 
the reduced density matrix vanishes. 

Let us first consider an initial state with no entangled 
qubit pairs, given as 



|^„(f = 0)) = 101111...). 



(13) 



Under the time evolution with the above Hamiltonian, 
the down spin located at the first site propagates, 
through spin-1 magnon excitations. Denoting a state 
with the down spin at site I by |^), the one-magnon 
eigenstates are given as \q) = J2i 4'i{q)\1)- Here the one- 
magnon eigenfunction for a chain of N qubits with open 
boundary conditions is 



lil) 



N + 1 



sm{ql),q- 



irn 
N+l 



,n=l,..7V, (14) 



and the one-magnon eigenvalues are e{q) = if cos q. The 
interaction term in the Hamiltonian does not figure in 
the magnon excitations in this sector with exactly one 
down spin. Now the time-evolved state at a time t can 
be written as 



(15) 



where the time-dependent wave function is given by 



11,' 



(16) 



The matrix elements of the reduced density matrix can be 
expressed in terms of the function ji^m- Let us consider 
the pair of qubits located at sites i and j We have B = 
\J^,l\^C=\J,.^\^,D = l~B-C,H = J,^,Jl,. 

Though we started with a state with no entanglement, 
in time pairs of qubits get entangled due to propagation 
of the down spin through magnon excitations. The differ- 
ent entanglement measures have been plotted in Fig. 1 as 
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FIG. 2: The concurrence and LE measures of entanglement 
are plotted against time for unentangled and entangled initial 
states for the last pair of qubits for an open chain of iV = 20 
qubits. 



functions of time for the pair of qubits at the end of the 
chain (i.e. i = N —1, j = N) foTc a. chain of N = 20 qubits 
with open boundary conditions. The entanglement mea- 
sures show time-dependent oscillatory behaviour, owing 
to the sinusoidal wave function. One interesting point is 
that the concurrence equals the upper bound on LE for 
one-magnon states. 

The impurity measure. My is largest for almost all 
time, as it has contributions from two types of entangle- 
ments, namely the entanglement of this pair with the rest 
of the chain and the mutual entanglement between the 
two qubits. The concurrence is for the most part smaller 
than the impurity measure, but it is always greater than 
or equal to the lower bound on LE. The entanglement 
measures for various nearest neighbour pairs of qubits 
are shown for two different times, Kt/U = 50, 150. All 
the measures show peaks for some pair of qubits, the 
peak shifting with the time. 
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It is important to note that all the measures have sim- 
ilar overall structures, for instance when the impurity of 
the pair of qubits is a local maximum, it is also the time 
when the pair is maximally entangled. Thus although 
purity (or lack of it) is not an entanglement measure, 
even in simple models as these there is a strong correla- 
tion between them. The LE as well as the concurrence 
share this property and hence we may say that usable en- 
tanglement is created with time at the end of the chain, 
even as the qubits become more mixed. 

Let us now turn our attention to an initial state with 
the first pair of qubits maximally entangled, and the rest 
of qubits in state |1), given as 

|Ve(t = 0)) = -i=|(01 + 10)11111). (17) 

The above state still has a definite number of down spins. 
Analogous to the unentangled state, we can write the 
state at a later time in this case as 

\Mt)) = ^Yl(7i,i+ 71,2)11)- (18) 

Now the nonzero elements of the reduced density matrix 
are given hyD = l — B — C, and 

Bij= ^|7^,l+7i,2|^ (19) 

Cij= il7,^l+7,,2|^ (20) 

H^J= |(7^+7;,2)(7.,i+7..2). (21) 

The envelope of the entanglement measures for the last 
pair of qubits as functions of time show similar behaviour 
as that of the imentangled state \ipu) shown in Fig. 1. We 
plot the different measures for the two cases of unentan- 
gled and entangled initial states in Fig. 2. The concur- 
rence of the last two qubits for the initial state \tpe) is 
smoother than the case when the initial state is is to- 
tally unentangled, namely \4'u), in which case the entan- 
glement shows rapid oscillations that repeatedly nearly 
vanish. The lower bound on LE, (LBLE) shows up the 
difference between the two cases, even in the envelope. 
A part of the entanglement in LE for the last pair of 
qubits shown in the figure can be ascribed to the entan- 
glement generated through dynamics, as in the case of 
the unentangled state, but a substantial part comes from 
the transport of initial entanglement of the first pair of 
qubits, through the Heisenberg dynamics. 

In Fig. 3 we show a density plot of (the lower bound 
of) the localizable entanglement between neighbouring 
qubits,I/iSi,;+i(f), as a function of I, the location of the 
pair and the time, for a chain of N = 20 qubits for an 
initial state with the first pair of qubits in a maximally- 
entangled state, that is \tpe)- The brighter regions in the 
plot correspond to larger entanglement for the given pair 
at the given time. We can see the front of bright "entan- 
glement" patch moving across the chain and reach the 



end linearly in time, getting reflected, and again moving 
back and forth. After one traversal th(^ bright patch is 
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FIG. 3: LBLEi^i+i{t) is plotted as function of / (along x-axis) 
and the time t, with brighter regions corresponding to larger 
entanglement, for a chain of A'^ = 20 qubits for the initial 
state IV'e)- The bright patch travels across the chain to reach 
the end. 



not so clearly seen, as entanglement begins to be shared 
by more than one neighbouring pair of qubits, simultane- 
ously. For the first passage across the chain, the location 
of the patch at time t is, I = vt, where v = K is the veloc- 
ity of the front. For a nonzero entanglement at a distance 
/, the pair of qubits located at 1,1 + 1 have to first be- 
come a part of a mixed state. This happens through the 
transport of the down spin located at the first site via 
one-magnon excitations. The time scale associated with 
the magnons is t = 1/K, and v = 1/t. 

In conclusion, we have investigated the dynamics of 

Heisenberg- XY spin chain using the various pairwise en- 
tanglement measures. Initial states with no entangled 
pairs of qubits, and with exactly one pair of entangled 
qubits are studied through the time-evolution, and the 
their dynamics of entanglement is juxtaposed. We have 
shown that the localizable entanglement bound captures 
the essential features of entanglement transport, and con- 
trasts better the dynamics of a state with no initially- 
entangled pair of qubits from the dynamics of a state 
with a pair of qubits maximally entangled. The entan- 
glement transport occurs through magnon excitations, 
and the first passage to the end of the chain of the pair 
entanglement front occurs with a uniform rate. 
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